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Abstract 

We report X-ray observations of the galai^ M87 and of a region In the 
Virgo cluster displaced 7* from the center of M87. We obtained X-ray spectra 
at these two locations with the Solid State Spectrometer onboard the Einstein 
Observatory. Emission lines were observed in both locations, indicating the 
presence of heavy elements at abundances ■> solar (to within a factor of 2). 
There is no strong abundance gradient, within our errors. We have detected a 
temperature gradient: T increases from 1.^/ keV at the position of M87 to 
T - 3.35 keV 7' way. There is lower temperature thermal emission at the 
center of M87 with T 0.6 keV, consistent with models for cooling flows in 
this cluster. We obtain a mass Infall rate of M ~ 20 Mq yr”^ . In addition to 
the thermal emission, we have detected a power-law component In the spectrum 
of M87, consistent with that observed by HEAO-1, Indicating that this com- 
ponent probably originates In the galaxy Itself. 

Comparison of our results with those of other observers and with 
theoretical models for this source indicates the presence of Intracluster gac 
having density 10“^ cm“^ and temperature 3 x 10^ K. 

Subject headings : galaxies: clusters of — galaxies: Individual (M87) 

— galaxies: Intergalactlc medium — X-ray: sources 
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I . INTRODUCTION 

The giant galai^ M87 In the Virgo cluster has been kncwn to be an X-ray 
source for many years. Recently It has become clear that the source is com- 
plex both spatially (Lawrence, 1978; Lea et al. , 1979; Forman et al. , 1979) 
and spectrally (Lea e£ al. , 1981) . In addition, attention has been focussed 
on this galaxy because of observations which may Indicate the presence of a 
black hole In the nucleus (e.g. Young et al. , 1978). In i. previous paper (Lea 
et al. , 1981) we presented spectral observations of this object using the A-2 
and A-4 detectors on board the HEAO-1 satellite. Ihese data Indicated the 
presence of a power law (or very high temperature thermal) component In addi- 
tion to a thermal component of kT = 2keV. We were unable to determine the 
spatial location of the spectral components we observed, due to thct large 
field of view of the detectors. However, the possibility of variability of 
the hard component led us to suggest that the nucleus of M87 was the location 
of the power-law component. 

The solid state spectrometer (SSS) onboard the Einstein Observatory has 
a 6* diameter circular beam, and offers the opportunity to obtain better spec- 
tral resolution, combined with spatial resolution, and thus to disentangle the 
various components of the source. In this paper we report SSS observations at 
two positions, centered on M87 and 7' offset from the center of M87. Observa- 
tions more distant from M87 proved Impractical due to the rapidly decreasing 
surface brightness. The distance to MS7 is taken to be 20 Mpc throughout, so 
that one arcmlnute is approximately equivalent to 6 kpc. 

II. OBSERVATIONS 

The observation of the central region of the galaxy M87 was obtained by 
the SSS in 1978 Jane 22-26, during which time a total of 16,500 seconds of 
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data were accumulated. The observed spectrum is complex (see figure 1) and no 
one single component model fits the data. Our basic model Is that of thermal 
emission from a hot, thin plasma (Raymond and Smith, 1977, RS), Including line 
emission from heavy elements. The Importance of the lines Is shown In figure 
2, which shows the ratio of lines to the "continuum" . 

Our analysis procedure was to attempt to fit a simple model, and to 
gradually increase ^the complexity until an acceptable fit was obtained. The 
amount of Ice accumulated on the detector was an additional free parameter, 
which was allowed to vary over reasonable values, as discussed by Holt et al. 
(1979). The results obtained were consistent with the expected value based on 
an empirical model for the rate of ice accumulation.' 

A two component RS model with the abundances for elements with z > 5 all 
allowed to vary together relative to H, gave a ot 221 for 58 degrees of 
freedom and was, therefore, unacceptable. Introducing a third thermal com- 
ponent, we found that the third temperature was so high as to be indistin- 
guishable from a power-law mimicking the free-free Gaunt factor, as indicated 
by our observations at higher energies (Lea et al. , 1981). For models 
including a power law of form F « E“® keV (cm^ sec keV)”^ , a is constrained to 
lie in the range 0.4 < a < 1.0. Adding a fourth component does not improve the 
fit. All the residuals appear to lie in the iron L band at around 1 keV. 

The best-fit model has a of 129, " 0.6 keV, Tjj » 1.4 keV, 

abundance of heavy elements (In solar proportions) approximately solar and an 
associated absorbing column 5 x lO^® < % < 2 x 10^^ cm"^ , plus a power law of 
energy index a ~ I with an associated absorbing column % '' 10^2 cm“^ . The 21 
cm hydrogen column in the direction of M87 is ~ (2 - 3) x 10^0 (Heiles 1976, 
Helles and Cleary, 1978). The cool component that we have observed itself 
contributes to the absorbing column, since about 15% of the oxygen at this 
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temperature has an Inner shell electron. We find % 3 x lO^O c«“2 

due to this material, where R Is the size of the cool component, defined below. 
Therefore, we do not consider our measured column for the thermal component to 
be In disagreement with the 21 cm value. However the measured value for the 
power law component Is SO times the 21 cm value and Is clearly discrepant. 

A formally acceptable fit (x^ * 87 for 55 degrees of freedom) was obtained 
by allowing the abundances to vary Independently. All were well within a 
factor of two of solar, except for Neon, which, with a value of about 6, more 
likely represents an Inadequacy In the atomic physics cot^iutatlons for the Iron 
L shell blend (Raymond, private communication, see also Smith, Mushotzky and 
Serlemltsos, 1979) than a physically significant abundance anomaly. It appears 
that the atomic physics calculations we are using (Raymond and Smith, 1977) do 
not accurately reproduce the detailed shape of the Iron L shell line group at 
'' 1 keV* Ne also has a transition in this region, and the shape of the emis- 
sion in this region can be changed by changing the Ne abundance. For this 
reason, the values in table I are for the case where the heavy element abun- 
dances are fixed ±o solar proportions. 

We have also observed a region In the cluster 7' from M87, during the same 
time period. We obtained 24,300 seconds of data at this location, R.A. 
12*’27®50®, declination 12°42'22’’. A simple thermal bremsst rah lung fit with no 
lines gave a 107 for 62 degrees of freedom. A better fit (x^ ■ 70 for 

61 degrees of freedom) was obtained by introducing heavy elements, in solar 
proportions, and the best fit for the ratio of heavy elements to hydrogen was 
0.95 solar. The Importance of line emission in this spectrum Is shown in 
figure 2. The gas tenq>erature for this fit was 3.35 keV, with 2.8 < T < 4.5 
keV at 90Z confidence, for 3 significant parameters (the temperature, abundance 
and emission integral). Ihe range of abundances determined In the same way Is 
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0.30 < X < 2.00. The amount of ice on the detector vae fixed at the value 
found for the M87 observation, although at this higher temperature (Hir results 
are not as sensitive to the amount of Ice. Introducing further components, or 
allowing the elemental abundances to change Independently did not produce any 
significant Improvement in the fit. 

III. DISCUSSION 

a) Abundance and Temperature Gradients 

The spectrum taken 7’ from M87 shows that X-ray lines due to heavy ele- 
ments are present at this location. The abundances at 7* from the center of 
M87 are consistent with those found at the center. Implying that no strong 
abundance gradients (i.e. greater than a factor of 2) are present. In con- 
trast, we find that the gas temperature is not consistent in the two observa- 
tions; the temperature is definitely lower in the observation centered on M87. 
This result is consistent with models for the X-ray halo which involve cooling 
in the inner regions of M87. The temperature which best fit the HEAO-l data 
(T 2 keV) is bracketed by the temperatures in the inner region (T ■ 1.4 keV) 
and farther out (T ■ 3.35 keV) . The lowest temperature (T ■ 0.6 keV) is repre- 
sentative of a very small amount of gas. The emission measure for this com- 
ponent is 4.6 X 10^^ cm~3 . Por a spherical volume of radius Ri^pc, this Implies 
a density of <n^ >1 /2 ~ o.6 ^ g mass of 6 x 10^ *^kpc^ 

pressure ^ 10“^ *^kpc~^ dynes cm""2 sufficient to confine the radio jet if Rj^pc 
< 0.4 (Owen et al. , I960). The temperature gradient in the Virgo cluster has 
also been estimated from the IPC image of M87 by Fabrlcant et al. (1980). Our 
measured temperatures in the inner 6' and 7* from the center of M87 are consis- 
tent with theirs. 

We have calculated the mass infall rate implied by our data as follows. 
Using the cooling rates in Raymond, Cox and Smith (1976), corrected to an iron 
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abundance of 3.2 x 10~^ , the cooling time at the S x 10^ K teoqperature we mea-' 
sure la 2 x 10® n~^ yr, where n is the gas density. Using our value of n 
given above» we find tcool ~ 3.3 x KP Rkpe^ 7^®* * total mass of gas 

6 X 10^ Rkpe^'^^ which implies M - 18 J^/yr, Independent of the else of the 
cool region. This result Is also Independent of the degree of cl\implng of the 
gas. We can also calculate the infall rate from our observations of the hotter 
component in the central region of M87, T - 1.4 keV ■ 1.5 x 10^ K. For this 
component the emission measure Is 2 x 10®® ■ <n^>V. Assuming that the source 
fills the SSS beam, we obtain >1 /2 « cm“® , t^ool ~ ^ ^ 

M ~ 9 x 10^ and M 23 yr"^ » In good agreement with our estimate from the 

e 

low temperature component. The computed value of M depends on the size of the 
1.5 X 10^ K region If it Is larger than the beam size, trlth M scaling approxi- 
mately as (1/23 kpc)~^ where I is the length along the line of sight of the 
emitting region. Ihls estimate Is sensitive to the iron abundance, and to the 
detailed physics of the cooling process. In particular. If the cooling Is 
Isobaric the computed M Is a factor of 2 lower. 

We have compared our results with the radiative accretion flow model of 
Mathews and Bregman (1978). Our measured temperatures of 1.4 keV and 3.3 keV 
over the central 6' and at 7' radius are reasonably consistent with their pre- 
ferred 30 H0 yr~^ accretion model, although our temperature gradient Is somewhat 
steeper (see figure 3). We have calculated emission measures for the model by 
approximating their density distribution by the formula n ■ 0.56 r“^ •** where r 
Is the radial distance In kpc. This formula is a good approximation to the 
numerical results for 1 kpc < r < 100 kpc. Our observed emission measure for 
the 1.4 keV gas Is In reasonable agreement with that predicted by the model 
(corrected for our aperture). The emission measure of 0.6 keV gas agrees very 
well with the model. On the other hand the emission measure predicted for our 
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bead at the 7' radius location is deficient by a factor of 2. Me conclude that 
there nust be more material at T 3 keV at larger radius than is predicted by 
the model (see also §b below). We note that this model implies that the tempera- 
ture should continue to Increaie outward. The emission measure of 8 keV gas 
in the model is 10®^ cm"^ , barely consistent with our previous limits from 
HEAO-1 (Lea et al. , 1981), although it is not clear how valid the model is at 
this large radius.' 

The radiative accretion model of Binney and Cowle (1981) provides a 
slightly worse fit, since the lowest observed temperature of S x 10^ K is never 
achieved in this model. The model of Takahara and Takahara (1981) is too hot 
everywhere, and does not fit our data. The temperature profiles for the three 
models, together with our observed tenq>eratures , are shewn in figure 3. We 
caution, however, that these models were designed to match to a hot Iptracluster 
gas (T ■ 10® K for MB and BC, T » 5 x 10^ K for TT). We have been unable to find 
any evidence for gas at 10® K in our data, (see also Lea e£ al. 1981 and §b 
below) . 

We can also compare our results with those of Canlzares et al. (1979, 1982) 
using the Focal Plane Crystal Spectrometer onboard the Einstein observatory. 
Their aperture is larger than ours, but may intercept about the same fraction of 
the emission measure of the 1.4 keV con^onent (based on the Mathews and Bregman 
model), because of its different shape (3' x 30' rectangular versus 6' diameter 
circular). They calculate an emission measure of about (3-4.5) x 10®** at 1.5 
X 10^ K and about 1.5 x 10®‘* at 5 x 10^ K. Their model assumes that there is a 
continuum of temperatures, and the curve of emission measure versus T is fairly 
flat at 10^ K. To obtain a better comparison with our simple, two component 
model we have integrated over the bandwidth 7.1 < log T < 7.3 to find the total 
emission measure in the peak of their distribution, and we find a value of 1.5 x 
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10^^ cm'*^ , which should be compared with our value at a coiq>arable tanq>erature 
(1.35 keV), 2 x 10®^ cm”^ . A somewhat smaller correction factor shmild be 
applied at 5 X 10^ K, which also brings the numbers from the t%iD experiments 
closer together. 

b) Inttacluster gas 

Fabrlcant e^ al. (1980) derived a density distribution from their surface 
brightness profile. The diffuse cluster background was subtracted from their 
data, so that this profile represents gas associated with M87. We have calcu- 
lated the emission measure expected In our 6* diameter beam from this density 
distribution.^ In the beam centered on M87 we find a predicted emission mea- 

^We have corrected the published profile for a typographical error, the cor- 
rect value of c' is 3.8 X 10"^. We thank Dan Fabrlcant for bringing this to 
our attention. We have also renormalized the density to account for the dif- 
ferent distance to M87 assumed by Fabrlcant et al . 

sure of ~ 3.3 X 10^^ cm"^ , as compared with our observed value (for both tem- 
perature components) of 2.5 x 10^^ cm”^ * Since the Fabrlcant et al. profile 
was derived assuming constant temperature, it overestimates the density In the 
cool. Inner regions. Therefore, these values are In reasonable agreement. At 
our second observed position 7' from M87 we infer from the IPC profile an emis- 
sion measure 6 x 10^** cm“^ , compared with the observed value of 7 x IC^** 
cm”^ • The total emission measure in the Fabrlcant e^ al. (1980) profile out to 
r - 300* (the OSO-8 field of view) Is 2.3 x 10^^ cm"^ , as compared with 6 x 
10^^ cm“3 observed by OSO-8 (Mushotzky et al. 1978). These results Imply that 
there Is additional gas with an emission measure of about 4 x 10^® cm"^ . This 
Is presumably the true Intracluster (as opposed to M87-assoclated) gas. If 
this gas gives rise to the large source observed by Lawrence (1978) and Lea el^ 
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al» (1979), it has a core radius of about 1” • 

We may derive a model for this Intracluster gas by picking a density dls* 

trlbution n - np(l + r^/a^)"^''2 where a Is the core radius, and choosing a and 

so that the total emission measure within 5" Is 4 x 10^^ cm~^ (to achieve 

consistency with the OSO-8 result) and the emission measure in the SSS beam at r 
« 7' la 1 X 10^** cm“^ (the maximum allowed for consistency between the IPC and 
SSS results). The result is a ■ 67.5', n^ - 9.6 x lO"** cm”^ ^ with the relation- 
ship ■ 9.6 X 10”** (a/l.l*)“^ , This size Is consistent with the HRAO-1 mea- 
surements. A larger size (and hence lower i^) Is excluded since both fields of 
view of the HEAO-1 A-2 experiment measured the same flux as was measured with 
the larger field of OSO-8 (Lea a^* 1981). A smaller size would lead to too 
much cluster emission in the Inner regions. The model Implies a higher surface 
brightness for the cluster component than that reported by Fabrlcant et al., 
although this surface brightness Is essentially uniform over the IPC field. 

Since the galactic background In this region Is very non-uniform (Lea et al. 
1979), making It difficult for the IPC to obtain an accurate background measure- 
ment for a source which fills Its beam, we do not consider this a critical con- 
straint. The cluster density of - 3 x 10“^ cm"^ , predicted using this model, at 
the position of M86 (1.25® away from M87) Is slightly less than the value 6 x 
10”** cm“^ estimated by Forman e£ al . (1979) at the position of M86. These 
results could be brought Into better agreement if the cluster center were 
slightly displaced towards M86, as may be Indicated by the optical data. 

We note that for this model, most of the HEAO-1 and OSO-8 flux is contri- 
buted by the "cluster” gas rather than "M87" gas. Thus the approximately 
solar abundance and 2 keV temperature measured by these Instruments (op . clt ) 
must apply to the cluster coo^onent. The temperature certainly cannot exceed 
the T ~ 3 keV value we measured at r - 7' . The density distribution resulting 
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from this analysis Is shovn In figure 4* While this model Is certainly not 
unique. It does Illustrate the Importance of the cluster component In 
Interpreting the X-ray observations of this region. The 3 keV temperature la 
consistent with the T-4V and T-emlsslon measure correlations derived from OSO-8 
data (Mushotzky e^ al. 1978) . Ihe total amount of Intracluster gas (In this 
model) Is 2.3 x lO^^ (a/1.1*)^'^^ Mq within 60', as compared with 1.9 x 10 ^^ 

In the Fabrlcant et al. profile, or 2 x 10^^ Mq (compared with 7.4 x 10^^ Mq) 
within 300*. The cooling time for this gas. If T ■ 3 keV, Is ~ 4 x 10^^ years, 
and hence this gas cannot participate In a cooling flow. 

c) The non-thermal component 

The large absorbing column associated with the power law component In our 
best-fit spectrum indicates that the most probable location for this component 
Is the galactic nucleus. This conjecture can be checked for consistency by com- 
paring the Intensity of the power law component observed In the two SSS fields 
with that observed by HEAO-1. In order to facilitate comparison of data sets, 
we fixed the power law Index at a * 0.73, a value typical of active galaxies 
(Mushotzky et al. 1980) and close to the predicted Inverse Compton spectrum, 
a ■ 0.79. For this model the normalization, N, of the power law was 0.013 ph 
(cm^ sec keV)“^ for the HEAO-1 data set taken in 1978 June, and 0.017 In 1978 
December. The SSS data from the central region were best fit with N « 0.012 
(see Table I) but could be fit by a value as low as 0.005, or as high as N •• 
0.018 In extreme models. Thus we conclude that it Is likely that the power law 
component originates within 3* of the center of Hc7. This result would be in 
accord with our Inference from the possible variability of this conq>onent and 
with the large observed value of %• If the temporal variation Is real, then 
c4t << 1'. Further, our observations In the cluster Itself 7' from the center 
of M87 do not require any power law component. This spectrum Implies a 1-60 keV 
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luffli'noslty of 9 X 10^^ erg . Using the luminosity function of Plcclnotti st 
si. (1982) we would have expected 1 source of luminosity > 9 x l(t*2 erg s~^ In 
the Virgo cluster. 

In our previous paper (Lea et al. , 1981) we also discussed this component 
In the context of a model for the hord X^ray emission of Inverse Compton 
scattering off the 3K background. We showed In that paper that the expected 
size of the Inverse Compton source would be about 16* x 12', the measured size 
of the Virgo A radio source at v > 80MHz. If this estimate Is correct, then the 
SSS should have detected only 15% of the Inverse Compton source. We conclude 
that at most thirty five percent of the hard X-ray emission observed by A-4 can 
be due to the Inverse Compton process, the rest being due to the nuclear source. 
This result has the effect of strengthening our previous limits to the field 
strength. If the source Intensity at the time of our observation were as 
observed In June 1978 by the A-4 experiment, ve find a stronger limit B > 2 x 
10“^ Gauss. This result is consistent with the value B - 2.5 x 10“^ Gauss 
estimated by Dennison (1980) from radio Faraday rotation data. However, If the 
source were In a high state, for example N <■ 0.0174 as we observed In December 
1978, then our limits to the B field are weaker, B > 7 x 10"^ Gauss. 

Einstein Observatory HRI measurements of M87 (Schreler 1981) Indicate that 
both the nucleus and the jet are sources of X-ray emission. If the possible 
variability we have previously noted (Lea et al. , 1981) Is real, then the 
nucleus itself must be the source of all or most of the power law emission. The 
observed column density In front of our power-law component (Njj ■ 2.4 x 10^^) 
Implies that the HRI source should be strongly absorbed. This column contri- 
butes an optical depth of almost 5 at 1 keV, sufficient to reduce the observed 
intensity by a factor of about 100. Our spectrum, with the best fit ■ 

0.012 E~‘ keV cm"^ sec"^ keV“^ and % • 2.4 x 10^2 contributes an average 
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flux of ~ 7 X 10"^ Jansky over the HRI bandpass, as compared with the observed 
flux of 6 X 10“' Jansky (Schreler 1981). These numbers are In reasonable agree- 
ment given the uncertainties In the calculation, but our source may also Include 
emission from the Jet, which la also within our beam. 

Possible models for such a nuclear source Include a Synchrotron Self 
Compton model (Maracher et al. , 1979; Lea et al» , 1981). In the context of this 
model, variations in the X-ray should be reflected in variations in the radio. 
There has been only one report of radio variability in the nuclear sourct, 

(Graham, 1971) which has appeared to be constant throughout all subsequent 
observations. This behavior Is puzxllng. Further searches for variability in 
both the X-ray and radio would be highly desirable. 

A variable X-rav source could also be produced In an accretion disc around 
a black hole. F.vidence has been produced for the existence of a black hole of 
mass ~ S X In the M87 nucleus (Young £t , 1878). This Interpretation 

of the data has been questioned, however, by Duncan and Wheeler (1980) and 
Dressier (1980) among others. The Eddington luminosity for such a black hole Is 
- b.S X 10** ’ ergs s“l, requiring an accretion rate of 10fl(Q/yr. We have 

seen that our observations of the thernwl sovirce are consistent with an accre- 
tion rate of 20 f\>/yr, which would give rise to an accretion luminosity of 
' 10“*’ ergs s“^, far In excess of what we see. nuis less than 10”*^ of the 
Inflowing material can end up In such a black hole, or else the accretion lumln- 
osltv does not emerge as X-rays. Tlie luminosity In othvT energy bands la also 
well constrained, e.g. L < A x 1(7*^ ergs/sec In the vlslhle (Dressier 1980), 
except In the EUV, where no limits exist. Much of the luminosity Is expected to 
be In this region If the disc temperature Is - 10^*’ - 1(7’ (e.g. Shakura and 
Sunyaev 1971, Pringle and Rees 1972, Novikov and Thorne 1973), but models with 
either a hot inner disc or disc corona predict l.^ ' Nmiv» then the predicted 
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exceeds our measured value* Note also that at 20 yr~^ the 5 x 10^ 
black hole would be formed In only 2.5 x 10® yr. 

There Is an Important caveat to the above discussion. The models we have 
discussed all assume disc accretion onto the black hole. If the accretion Is 
essentially spherical (l.e. if the accreted material has very low angular 
momentum) then the efficiency for producing energy from the accretion flow may 
be much reduced. For models which Include heating effects due to turbulence 
and/or magnetic field effects (Ipser and Price, 1977, Maraschl e^ al. 1979) 
high efficiencies are regained. In particular, Maraschl et al. calculate the 
luminosity from a 10® black hole accreting at few x 10 1% /yr and find 

total luminosities ~ few x lOf*® ergs s"^ , with a large fraction of this In the 
X-ray band. These models exhibit a wide range of values for the fraction of 
total luminosity emitted in the X-ray range, but all violate either our 
measurements or Dressier' s (1980) optical measurements. 

We conclude that if a massive (5 x 10® Mq) black hole exists at the center 
of M87, the infalling gas mast be prevented from reaching It, unless the accre- 
tion flow Is essentially spherical and non-dlsslpatlve. 

d) The fate of the accreted material. 

If the accreted matter does not reach the black hole, where does It go? 

Our observed absorbing column of % '' 10^2 cm“^ for the power law component 
implies the existence of 320 r 2 pc of cold material In the nuclear regions of 
M87. This amount Is negligible for any reasonable values of the radius R of the 
nuclear region. Dressier (1980) has concluded that the central luminosity spike 
In M87 Is a cluster of stars. However, he estimates the age of the cluster to 
be > 10® years. If we take a conservative value of M « 10 yr“^ then, 10^® 
of gas would have accumulated In the core since those stars formed, and this gas 
is not observed in either neutral or ionized form. Even If we speculate that 
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the gas forms Into low mss stars only, so as to be consistent with Dressier* s 
observed spectrum In the central spike, then apart from the unusual Initial mss 
function (IMF) this Implies, we also require a mss to ll^t ratio 1000 for 
these stars In order to have a star cluster of total luminosity 1(]^ • This 

seem unlikely. 

lUE observations of M87 (Bartola e£ al. 1980, Perola and Tarenghl 1980) do 
Indicate the presence of som blue stars In M87. These my be either horizon- 
tal branch stars or young blue stars. Even if they are young blue stars, our 
comparison of the M87 flux with that from local stars (Brune, Mount and Feldmn 
1979) Indicates that at most 10^ could be present in stars of mss 10-30 . 

Assuming a Salpeter IMF, this Implies 10^ of stars In the mss range 0.1- 
30 . These stars must have formed over the ~ 10^ yr lifetime of a 10 

star, Implying a gas consumption of about 1 yr"^ , or less than one tenth of 
the total mss Inflow. A steeper IMF of the form N(M) « could accomodate 
the total Inflow rate at 10 yr“^ in star formtlon. Note, however, that 
this implies that all of the observed UV light Is due to young stars of 10-30 
M0 and none is due to blue horizontal branch stars, which Is perhaps unlikely. 
For a 20 yr”^ Infall rate, a steeper IMF would be required. In this 
context, we should mention that there have been suggestions that only low mss 
stars form in elliptical galaxies (Jura, 1977) or in a cooling flow (Sarazln, 
private communication) . If this Is the case, the UV observations are much less 
constraining and all of the infalling gas my form into stars. 

Recent work by Gunn, Stryker and Tinsley (1981) also Indicates that star 
formtlon Is occurring in giant elliptical galaxies. For M87, they find that 
the observed blue light may be due to "young stars", implying a star formtlon 
rate of 1-2 M0 yr”^ , or the blue light my come from evolved "blue straggler" 
stars, implying that 10-30% of the stars in the galaxy are of this type. In 
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the 'latter case, we estimate a rather uncertain star formation rate of about 20 
Me yr"^ . These authors favored a flatter (rather than steeper) IMF, N(M) « 
with a < 1, and an upper mass limit of 2 in the stars formed. 

The radio source Virgo A and the observed X-ray /optical/radio jet imply 
that ~ 5 X 10^^ ergs yr“^ of energy are being ejected from the nucleus of M87. 
Mathews and Bregman (1978) have noted that this power is comparable to the 
power in the accretion flow, and they speculate that the infalling material is 
re-ejected into the radio source, perhaps in the Jet. Recent work by DeYoung 
(1981) indicates that matter may be entrained in the jet and carried outward to 
the radio lobes. The mass is entrained in the nuclear regions where the inter- 
stellar density may be relatively high. The rate at which matter may be 
removed from the nucleus by this process is 

M - 60Rj 6 X n^ X Cg ^^/yr. 

where Rj is the Jet radius in kpc, r»o is the gas density in the entrainment 
region, 6 Is a parameter, 6 1, and eg is the sound speed in the entrainment 

region In units of 1000 km/sec. If entrainment occurs in the region of our 
lowest temperature con?>onent, we may calculate values for n^ and eg. DeYoung 
argues, however, that because o( Interaction with the jet the material in the 
entrainment region next to the jet may have cg ~ 1. The calculated value of n^ 
depends on the radius of the central region. For a radius of ~ 0.1 kpc, Oq 
20 cm~^ • The jet radius Is obtainable from the VLA map (Owen et al. 1980); we 
take Rj < 50 pc to find 

M < 36 cg /yr 

Tlius this mechanism is not sufficiently efficient to remove the accreted 
material unless the entrainment region lies Interior to our observed 0.5 keV 
gas and has n^, >> 20. Recent VLA observations of the jet (Owen e^ ajl 1980) 
Indicate a density of perhaps 0.1 cm~^ In knot A. Even If material of this 
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density were present throughout the jet» the maxltaum mass outflow rate would 
be 6(v/c) t^/yr where v Is the speed of the material In the jet, aich less 
than the required 10-20 Vq/yt for any reasonable values of v/c> Thus, If the 
DeYoung mechanism Is operative, the entrained material must remain unmlxed 
with the jet magnetic field. 

Another possibility Is of the existence of some sort of Invisible ejec- 
tion, either a continuous wind taking place concurrently with the Inf all and 
Interleaved with It, like segments of an orange, perhaps, or sporadic violent 
ejection. In this latter scenario, the most recent outburst must have oc- 
curred fairly recently (within about a million years) so that the amount of 
material accumulated since the outburst remain Invisible, but long enough ago 
that no obvious vestiges of Its occurrence remain. In this picture the large 
radio source Virgo A and the jet may be the long-term residuals of these ex- 
plosive outbursts. The ultimate cause of the outbursts Is a mystery. 

One further possibility is that our inflow calculations are in error. 

The inflow need not occur if there exists a heating mechanism to ba?.ance the 
cooling of the gas and prevent inflow. Such an equilibrium state Is not only 
improbable but may be unstable. 

e) The origin of the infalling material 

Previous models of this object have assumed that the material we see 
emitting X-rays has come from the member galaxies of the Virgo cluster. 
However, it is worth noting that the 10-20 yr“^ which we see accreting onto 
the M87 nucleus could be provided by M87 itself, if the galaxy is massive. If 
the mass loss rate from stars is 10“^^ - 10“^^ (^^ yr)”^ (e.g« Faulkner and 

Freeman, 1977), then we need a total stellar mass of 10^ ^ - 10^^ to provide 
the necessary material. There have been several estimates of the mass of M87 
in this range (e.g. Mathews, 1978, Fabrlcant et al. , 1980, but see also Blnney 
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and Cowie, 1981)* Of course most of this mass is invisible and we have no 
evidence to suggest that this mass is shedding gas« Certainly if low mass 
stars, neutron stars or black holes are responsible for the large galai^ mass 
estimates, no mass loss from these objects is expected* Faber and Gallagher 
(1976) have estimated the mass loss rate in terms of the galaxy luminosity. 

Using their value and the luminosity of M87 from Oemler (1976) we find a total 
mass loss rate of 3 yr~^ plus or minus a factor of 3. Given the 
uncertainties in our knowledge of all these numbers, it is certainly possible 
that M87 is simply eating itself. Sparke (1981) has recently discussed models 
of this type. If this is the correct explanation for M87, then similar objects 
may be found which are not members of rich clusters. Kriss et al. (1980) have 
reported that three D or cD galaxies in small groups are x-ray sources. Perhaps 
these are also self -cannibals. 


IV. CONCLUSIONS 

We have observed the X-ray spectrum of M87 and a nearby region in the 
Virgo cluster. From our spectra we conclude that the power law component 
previously observed in this source most likely originates in the galai^ M87 
itself. This result allows us to improve our previous limit to the magnetic 
field strength in Virgo A to B > 2 x 10"^ Gauss. The non-thermal X-ray source 
that we have observed is not unusual for an active galaxy, and probably does 
not result from some peculiarity in M87 Itself, but is an example of a more 
general phenomenon. 

For the thermal source we conclude: 

(i) We see no evidence for strong abundance gradients in our data. 

(il) We see a strong temperature gradient, T increasing from-' l.A keV 

in the inner 6* to T - 3.35 keV at r ^ 7' . In addition, there is a 
small emission measure component at T ~ 0.6 keV in our M87 data. We 
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stress that this is the simplest model that fits our data* Such a 
model Is probably Indicative of a continuum of temperatures ranging 
from 3*35 keV dovn to at least 0.6 keV, with emission rieasure 
decreasing as a function of temperature* Such models are discussed 
by Canlzares et_ al * (1982) . 

(Ill) Our results are consistent with previous observations using the IFC 
and with the radiative accretion model of Mathews and Bregman 
(1978), and Imply M " 10-20 yr“^ * Infall at this rate Is 
difficult to reconcile with the presence of a black hole of mass 
5 X 10^ in the nucleus of M87, and poses severe problems as to 
the ultimate disposition of this material unless the accretion is 
essentially spherical and with a very low (e < 10”^) efficiency of 
photon production. However, by pushing all the numbers to the limit 
we can dispose of the material by a combination of star formation 
and jet entrainment. In particular, we require that the stars 
formed be lew mass, (< TMq ) . 

(iv) Comparison of the SSS data, the IPC data (Fabrlcant et al. 1980), 
the HEAO-1 data (Lea et al. 1981) and the OSO-8 data (Mushotzky et 
al. 1978) implies the existence of true Intracluster gas having a 
core radius of about 1“ , a central density of 10“^ cm“^ and 
temperature < 3 keV. 

(v) The abundance ratios in this cluster are consistent with solar 
values. 

We acknowledge helpful discussions with C. Canlzares, who provided us 
with his results prior to publication. During the course of this work SML was 
supported by NSF grant AST 78-21070, NASA contract NAS 8-33357 and NASA grant 
NAG 8385. SML acknowledges the hospitality of the Aspen Center for Physics 
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TABLE 1 

SPECTRAL PARAMETERS FOR M87 AND THE VIRGO CLUSTER 



PARAMETER 

M875 

7' FROM M8T5 

Tl 

0.58 keV 

- 

EM(Tj)‘ 

X 10‘" 

4.3 X 10®'' 

< 7 X 10®^ cm"^ 

"2 

1.36 keV 

“ • io 

3.35 ± keV 

EM(T2)^ 

, +.48 ,.6 5 -3 (3 ) 

1.77 _ X 10 cm ' ' 

7.2 tj;” X c.-' 


Nn(thermal 

component) 

Abundances:^ 


2.0 


+.5 

-.3 


X 


10 


.6 5 


cm 


-3 (*♦) 


1.3 


+0.7 

- 0.8 


X 


10 ^^ 


-2 

cm 


1.1 


+1.0 

-0.4 


1.7 


+ 1.0 

- 0.6 


X 


10 ^' 


-2 

cm 


0.95 


+1.0 

- 0.6 


Power law normalizatlonC^ >'* ) 0.012 ± 0.002 

Nh (power law con 5 >onent)(‘* ) 2.4 ± 1 x 10^2 cm“2 


None required 


^ At solar abundance. 

2 Abundances are relative to solar values as given in Allen (1973), exept Fe, 
for which we take a solar value of 3.2 x 10”^ relative to H. 

^ Power law index free to vary. 

^ Power law index fixed at 0.75. 

^ All errors are 90% confidence for two significant parameters (x^ + 4.6). 
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Figure Captions 

Figure la: The pulse height data for M87 from the SSS for the central pointing 
position. The solid histogram Is the best fitting model to t. ' 
data. 

Figure lb: The photon spectrum of the central position. Sharp spectral 

features show the Instrumental resolution of the SSS (160eV FWHM). 

Figure 2: The ratio of the SSS pulse height to "continuum” as a function of 

energy for the central postlon (solid line) and the position 7* off 
center (the dotted line) . The expected energies of the most 
prominent lines are indicated. The "continuum" has been determined 
separately for the central and offset positions, and Is a smooth 
exponential which passes through the lowest points In the data. 

Figure 3; Observed temperatures for the M87 source are shown along with 

models by Mathews and Bregman (1978) (MB), Binney and Cowie (1981) 
(BC) and Takahara and Takahara (1981) (TT). 

Figure 4: One model of the gas density distribution In the Virgo cluster, as 

deduced from the x-ray data. The dashed line Indicates the total 
density, while the two solid lines Indicate the (somewhat 
artificial) division Into "M87" and "cluster" components. The bar 
and arrow at the top of the figure indicate the extent of the IPC 
observations of this source. The cluster component dominates 


outside of the IPC field. 
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